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ON  THE  USE  OF  REPLACEMENTS  TO  E}CTEND  SYSTEM  LIFE 


by 

C.  Derman,  Columbia  University 
G.  J.  Lleberman,  Stanford  University 
S.  M.  Ross,  University  of  California,  Berkeley 


1.  Introduction 

The  literature  dealing  with  the  scheduling  of  replacements  and/or 
Inspections  is  very  extensive.  See  [2]  for  a  comprehensive  set  of 
references .  . 

In- th±*  -paper  we  are-  concerned  with  the  following  question.  A  system 
has  one  vital  component  for  which  there  are  n  spares.  Whenever  the 

r  C.  *  £ 

vital  component  fails,  the  system  fails,  mr  would  like  to  schedule  the 
replacement  of  the  vital  component  with  the  spares  so  as  to  prolong  the 
life  of  the  system  as  much  as  possible.  This  problem  can  be  generalized 
to  where  there  are  several  components  in  the  system  and  the  scheduling  of 


replacements  refers  to  each  of  the  components.  deal  mainly  with  the 


first  question  but  treat,  to  some  extent,  a  special  case  of  the  second. 


\ 


In  the  next  to  final  section  we  consider^  generalization  which  allows  for 
the  vital  component  to  fail  a  fixed  number  of  times  before  causing  the 
system  to  fall. 


<  -- 


2.  One  Vital  Component  Problem 


Assume  a  one  component  system.  Let  F,  with  F  *  1-F,  denote  the 
life  distribution  of  the  component.  He  assume  that  its  mean, 

00 

H  -  J  F(t)dt  <  -  , 

0 

and  that  F  has  a  continuous  density  f  with  support  (0,T),  0  <  T  <  ®, 
and  we  let 


r(t)  -  f(t)/F(t) 

denote  its  failure  rate  function. 

Suppose  the  system  is  put  into  operation  at  t  ■  0,  using  an  initial 
component.  There  are  n  statistically  identical  spares  that  can  be  used 
sequentially  as  replacements.  The  system  fails  whenever  the  component 
that  is  currently  in  use  falls.  The  idea  is  to  judiciously  replace  the 
component  in  use  with  an  available  spare  in  order  to  extend  the  life  of 
the  system.  He  assume  that  once  a  component  has  been  removed  it  cannot  be 
used  again.  Let  S  denote  a  replacement  schedule,  L  the  life  of  the 
system,  and  EgL  the  expected  life  of  the  system  as  a  function  of  S. 

The  problem  of  concern,  here,  is  to  determine  the  schedule  S,  and  (or) 
some  of  its  properties,  that  maximizes  EgL. 

As  will  be  evident  there  always  exists  an  optimal  schedule,  although 
it  may  not  be  unique.  Let  v^  denote  the  expected  life  of  the  system 
given  an  optimal  schedule  of  replacement  of  n  spares.  If  n  ■  0,  there 
is  only  one  schedule;  viz,  the  one  that  never  replaces.  Therefore, 


v0  "  14 


If  n  1,  vq  satisfies  the  optimality  equation 


v  -  max  $  (x)  , 

n  0<x<T  *» 


where 


♦  (x)  -  /  tf(t)dt  +  F(x)(x  +  vn  .) 
n  Q  n-I 


I  F(t)dt  +  F(x)  v  . 
/»  n-i 


Since,  by  assumption,  ix  <  ®,$  (x)  is  continuous  at  all  x  £  [0,T]. 

n 

Thus,  (x)  can  always  be  maximized,  although  not  necessarily  uniquely, 
n 

at  x  -  x  (say).  The  value,  x  ,  is,  in  words,  the  optimal  length  of 
n  n 

time  to  use  a  component  before  replacing  it  when  n  spares  are 

available.  If  x  -  T,  then  the  component  is  never  replaced  and  v  ■  jx. 

n  n 

If  xq  ■  0,  a  component  is  Immediately  replaced,  in  idiich  case  vn  »  vn_i‘ 
It  should  be  kept  in  mind  that  the  system  may  fall  before  all  n  spaces 
are  used  even  though  a  schedule  is  optimal. 


Proposition  1 :  vn  is  nondecreasing  in  n. 


Proof:  v  >  ♦  (0)  ■  v  , ,  n  >  1. 
— *  n  — ■  n  n-1  ■ — 


Proposition  2:  If  v,  ■  v„,  then  v  ■  v„,  n  >  1. 
- -  l  u  n  u  ~ 


Proof :  From  (1)  and  (2)  we  obtain,  for  n  2, 


-3- 


Since  vn-1  -  vn_2  0,  by  Proposition  1  and  0  _<  P(xn)  <  1,  the 
proposition  follows. 

Since  ^(0)  »  <^(T)  ■  M->  and  from  (1)  and  (2),  it  follows  that 
v^  ■  Vq  if  and  only  if 

x 

(4)  p>  J  F(t)dt  +  F(x)p  *  0<x<T  , 

0 

or,  equivalently,  if  and  only  if, 

p  <  /  ^ ^  t  )dt  o  <  x  <  T 

“  x  F(x) 

The  above  condition  is  often  called  the  new  worse  than  used  in  expectation 
(NWUE)  for  it  states  that  the  expected  additional  life  of  a  used  vital 
component  of  age  x  is  at  least  as  large  as  the  expected  life  of  a  new 
one. 

By  considering  the  derivative  of  the  right  side  of  (4)  we  can  see 
that  if  f (0)  -  r(0)  <  — ,  then  v ^  >  vQ.  On  the  other  hand  if  F  is 
NWUE  then  x^  can  equal  T  and  v^  ■  v^.  We  shall  assume  for  the 
remainder,  unless  stated  otherwise,  that  F  is  not  NWUE  and  so  v.  >  v~. 


Proposition  3:  v  >v  , ,  v  -  v  ,<v  ,-v  and 

- L -  n  n-1  n  n-1  n-i  n-2’ 

0<x  <x  ,  <  T,  n  >  2. 
n  —  n-1  — 


Proof ;  Since  v^  >  v^,  x^  *  0  or  T.  In  a  manner  similar  to 
showing  (3)  we  have 


-  Vl  >  '<Vl«Vl  •  Vn-2>  •  "  >  2 


If  x  *  T  for  some  n,  then  v  *■  n  ■  v„  which  is  contrary  to 
n  n  0 

proposition  1  in  view  of  the  hypothesis  that  v  >  v  .  Therefore,  x  <  I 

10  n 

for  all  n  ^  1,  and,  hence,  F(xn)  >  0.  Thus,  it  follows  inductively  from 

(5)  that  v  >v  ,,n>l.  If  x  “0  for  some  n  >  1,  then  v  *•  v  , , 
n  n-1  —  n  —  n  n-1 

a  contradiction  of  vrtiat  has  Just  been  shown.  Thus,  x^  >  0,  F(xn)  <  1, 

n  ^  1.  From  (3)  and  the  positivity  of  v^  -  vr_^,  it  follows  that 

v  -v  ,<v  ,-v  _,  n  >  2.  Finally,  from  (3)  and  (5)  and  *ft\at  has 
n  n— 1  n-1  n — Z  — 

Just  been  shown,  we  have 


F(xn)  >  F(xn_x)  ,  n>2  , 

from  vdiich  x  <  x  , ,  n  >  2 
n  —  n-1  — 

As  previously  remarked,  x^  need  not  be  unique.  However,  let  {x^} 

denote  the  set  of  values  for  which  the  maximum  of  4*  is  obtained.  The 

n 

proof  of  Proposition  3  shows  the  following 


Proposition  A:  x^  satisfies 


r(V  *  ■  "  i  1 

n-1 


and  r  cannot  be  decreasing  at  x^. 


Proof:  From  proposition  3,  x  is  an  interior  point  of  the  interval 
-  n 

[0,T].  On  differentiating  we  get 


$’(x)  -  F(x)  -  f(x)  v  , 
n  n  1 

-  F(x)(l  -  r(x)  v  .)  . 

n  i 


Thus,  <t£(x)  ■  0  only  if  r(x)  ■  l/vn_^.  However,  if  r(x)  -  -  and  x 

n  n  vn-l 

is  a  point  of  decrease  of  r,  then  is  increasing  at  x  so  that  x 

cannot  be  a  relative  maximum. 


Corollary:  x^  <  xn-1»  n  >  1. 


Proof:  Proposition  3  asserts  that  v  .  >  v  „  and  allows  the 
-  n-1  n-2 

possibility  that  x  -  x  However, 

n  n— i 


n-1 


r(xn> 


•  r(Vl> 


-  l/v 


n-2 


is  impossible. 


-6- 


The  second  derivative  of  ♦"  is 

□ 

<T(x)  -  -f(x)  -  f * (x)  v  . 
n  ni 

Thus,  If  x  is  a  root  of  r(x)  -  — —  and  f'(x)  0,  then  x  Is  a 

n-1 

relative  maximum  of  $  • 

n 

If  F  Is  IFR,  then  x  Is  the  unique  root  of  r(x)  ■  1/v  , . 

n  n-l 

Because  in  this  case,  f(0)  ■  r(0)  <  1/n.  In  fact,  from  Proposition  3, 

f(0)  -  r(0)  <  1/v  ,  n  >  1. 

n 

If  r  is  unlmodal  and  has  two  roots  to  r(x)  -  1/v  . ,  then  x  is 

n— l  n 

the  smaller  of  the  roots;  if  r  has  the  so  called  "bathtub"  shape,  then 

x  is  the  larger  of  the  two  roots • 
n 

Proposition  5:  lim  v  ■  v  is  finite  or  infinite  according  as  f(0) 
n -*■<*>  n 

is  positive  or  zero.  Although  x  is  not  necessarily  unique, 

n 

11m  x  ■  x*  >  0,  is  a  unique  value. 
n-Mo  n  — 

Proof:  If  v,  -  v„,  then  v  *  p,  n  >  1,  by  Proposition  2.  Suppose, 
— _  l  u  n  •“ 

for  the  remainder  of  the  proof,  that  v^  -  Vq  >  0.  Since,  by  Proposition 

1  v  >  v  then  11m  v  exists  and  is  either  finite  or  infinite.  By 
n+1  —  n  n-*»  n 

Proposition  3  and  the  corollary  following  it,  it  is  easily  shown  that 

lim  x  exists  and  equals  a  unique  non-negative  value  x*  (say).  Suppose 
n-*»  n 

f(0)  ■  r(0)  >0,  x*  ■  0,  lim  v  -  For  e  >  0,  choose  6  >  0 

n-*°  n 

sufficiently  small  so  that  r(x)  ]>  e  if  x  <  6.  If  v^  >  1/e,  then 
1/v  <  e,  «v  n  >  which  implies  that  x>6vn^N,  a  contradiction 


that  x*  »  0.  Thus,  under  these  circumstances  lim  v  -  00  is 

n 

impossible.  Suppose  f(0)  ■  r(0)  >  0,  x*  >  0.  From  Proposition  2 


Thus 


vn  '  Vl  i  F(xn)(Vl  •  vn-2> 
<  F(x*) (vn-1  -  vn_2) 


<  Fn  1(x*)(v1  “  vQ)  ,  n  >  2  . 


lim  v  -  v.  <  l  Fn  1(x*)(v.  -  vQ) 
n  -+<x>  n  1  —  n«2  ^  ^ 


Ii**l  (V 

F(x*)  V 1 


V0> 


<  ®  • 


Suppose,  now,  that  f(0)  ■  0.  Consider  the  schedule  S  having  equal 
replacement  intervals  of  length  y  with  ny  *  ^F(~y)^  e(lual  to  the 

number  of  available  components,  including  the  initial  one.  Let  L  denote 
the  life  of  the  system  under  S  .  Then  it  is  easy  to  show  that 

ny 

lim  E(L)  -  -  . 

y+0 


Since  v  >  E(L)  and  v  increases  in  n  it  follows  that  lim  v  ■  ® 
ny  n  n-*»  n 

If  r(0)  >  1/n,  then  x*  >  0.  For  it  is  trivially  true  if  x^  ■  T; 


and  if  x^  <  T,  then 


r(x*)  ■  lim  r(x  ) 
n-*®  n 

-  lia  1/ v  . 

n>°°  ®" 1 

-  1/v 

<  i/*0 

■  1/p  . 

The  following  Is  a  sufficient  condition  for  x*  to  equal  0. 

Proposition  6:  If  F  is  a  strictly  new  better  than  used  (NBU(S)) 
distribution  in  the  sense  that 

^  <  F(t) ,  for  all  x  >  0  ,  t  >  0  , 

F(x) 

then  x*  ■  0. 

Proof:  Suppose  F  is  NBU(S)  and  x*  >  0.  For  each  n  let  S  be 

n 

the  schedule  (xn/2*  Xn/2*  *n-2*  ***»  !•••»  the  first  two  replacement 

intervals  of  are  of  length  \/2  and  optimal  thereafter.  Let 

Gn(t)  -  F(t)  ,  t  <  xq^2 

•  F<V2)m-,n/2>  ,  t  >  xn/2  . 

Note  that 

Gn(t)  >  F(t)  ,  t  >  xn/2  . 


•  ■1 

n 


• .  « 


Let  v 


E_  (L).  Ve  have 


V 

D 


°  G  (t)dt 
n 


+  G  (x  )v  . 
n  n  n-2 


Denoting  G  -  11m  G  and  letting  n  ■*  •»  we  have 
n+<»  n 


11m  v  -  v  *  /  (G(t)  -  F(t))dt  +  (G(x*)  -  F(x*))v 

n-»®  n  0 


>  0  , 


a  contradiction  that  11m  v  <  v. 

n-*x  n 


Remark:  The  condition  that  F  is  NBU(S)  is  weaker  than  that  F  is 
strictly  IFR,  i.e.,  that  r(t)  is  strictly  Increasing. 

Examples: 

(i)  Uniform  <0,1).  Here, 
f(t)  -  1 

F(t)  -  t  ,  for  0  <  t  <  1  , 

r(t)  -  1/1-t 

vQ  -  |i  -  1/2  ,  vx  -  5/8  . 


Since,  f(0)  >  1,  v  < 


From  Proposition  4 


n-i. ' 


‘  Vl>  +  Vl> 


■z -  +  V  , 

2  n-1 


-  1/2(1  +  v*  )  ,  n  >  1 

n  l  — 


Some  values  of  x  and  v  will  appear  in  Table  1.  Since  F  Is 
n  n 

NBU(S)  we  see  that  x*  •  0  and  v  -  ■  1* 


(ii)  Linear  (0,1);  specifically, 


f(t)  -  2t 


F(t)  -  t  ,  for  0  <  t  <  1 

/  %  2t 
r(t)  “  - 2 

1-t 


Vg  ■  u  “  2/3  . 


Since,  f (0)  -  0,  v  -  x*  •  0  and 


-  rr 

1-x 

n 


-  1/v 


n-1  » 


leading  to  the  quadratic  equation 


x  +  2v  .x  -  1  -  0  , 

n  n-1  n 


for  which  the  only  positive  root  is 


"  r™v  ‘  ^  .  *  M  »  >  "  ^  -  • 


-  *  ^  V—  ¥  -  * 


*'/• 


Thus, 


K  “  "V  ,  +  /^77  • 

n  n-i  n-i 


v  -  2  /  n  t2  dt  +  <1  -  xh(xn  +  v  .) 
n  q  n  n  n—  i 


'  \  +  -  ‘.XV  +  Vl>  ’ 


The  values  of  and  can  be  computed  recursively,  commencing 

with  Vq  -  2/3. 


(ill) 


f(t)  -  2(1  -  t) 


F(t)  -  2t  -  t  ,  for  0  <  t  <  1  , 


"I=t 


vQ  -  p  -  1/3  . 


J 

1 


Here,  f(0)  ■  2;  hence  v  <  ®.  F  is  IFR  with  r(0)  <  —  ;  thus 


a 

% 


*co 


n'  1-x 


uV. 

•v.\ 

i 

I 


i 


K.v.y: '-v 

fcv&Iv: 


1  •  e  • , 


Also, 


1/v 


n-1  , 


x  ■  1  -  2v  , 
n  n-1 


v„  -  x2(l  -  2/3  x  )  +  (x  -  i)2(x  +  vn  .)  . 

n  n  n  n  n  n— i 
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The  values  of  x  and  v  can  be  calculated  recursively.  To  find  v,  it 
n  n 

is  easier  to  express  v  in  terms  of  x*  and  solve  for  x*  in  the  latter 
equation  after  having  let  n  -►  ®.  Out  of  this  procedure  comes  the 
equation 


x*2( 1/2  -  x*/3)  -  0 


from  which 


x*  -  0 


is  the  only  acceptable  solution.  Thus, 

v  -  1/2  . 

(iv)  Two  exponentials  in  parallel. 

Let 


f(t) 

F(t) 

r(t) 


v 


0 


t  t 

2e"  (1  -  e"  ) 

—  ^  2 

(l-e  )  ,  0  <  t  <  » 

t  t 

*  fi-«~  ) 

t  7 

i-(i )2 

I*  -  3/2  . 


This  case  arises  when  the  component  comprising  the  system  is  a  module 
consisting  of  two  independent  components,  each  component  having  an 
exponential  distribution  with  mean  one. 


3.  Equal  Interval  Replacement  Schedules 

In  this  section  we  restrict  the  class  of  schedules  to  those  having 

replacement  intervals  of  equal  length.  Let  4»  (y)  denote  the  expected 

n 

life  of  the  system  when  y  is  the  length  of  the  interval.  Let 

4,  m  max  (p  (y) 
n  y  n 

Let  y  denote  a  value  of  y  such  that 
Jn  J 


It  is  seen  from  (7)  below  that  such  a  value  exists  and  is  positive. 
Notation  from  Sections  1  and  2  will  be  continued. 

First,  we  note  the  recursive  relation 


4- (y)  •  f  tf(t)  +  F(y)(y  +  4>  ,( y>) 

n  0  n  1 

-  /  F(t)dt  +  F(y)  4'  ,(y)  ,  y  <  0  . 

a  Ha 


Also,  by  considering  the  difference  between  having  n  and  n-1 
spares,  we  have 


y 

4»  (y)  -  *  i<y)  “  yn_1<y>{/  tf(t)dt  +  F(y)(y  +  -  p} 

n  n-1  0 

y 

-  Fn_ld(y){J  F(t)dt  “  F(y)n} 

0 

■  Fn”1(y){^1(y)  “  p}  . 


(6) 


On  summing  (6)  we  get 


(7)  4>n(y)  -  4>0(y)  -  <J<n(y)  -  p 

y 

-  (/  F(t)dt  ~  F(y)p)(- 
0 

■  (*»>  - 

Proposition  7:  If  >  vQ,  then  4>q  >  <j»  , 

v  -  vQ  -  |i  then  4«n  ■  for  n  I. 

Proof :  From  (7), 

♦«  *  “  i  W  '  “ 

*  nxj  (vi  ‘  “> 

>  0  • 

i 

! 

That  Is,  <1»  (y  )  >  p,  n  >  I,  which  from  (7)  implies 
n  n  — 

♦l<yn)  >  P  .  n  >  1  . 

Thus,  from  (6) 

♦»  '  Vi  i  VVi>  •  Vi'Vi’ 

■  - 

>  0  . 


_  1  c.  _ 


Suppose  v^  ■  p.  The  above  argument  shows  that  <Ji  2  n.  However,  since 

v^  2,  and  vq  ■  p,  n  >  1,  by  proposition  2,  it  follows  that 

4»  -  P,  n  >  1. 

n  — 

A  reasonable  conjecture  is  that  for  any  set  of  maximizing  values 

{yn>,  yQ+1  <  yQ,  n  2  1«  We  prove  that  this  is,  in  fact,  true. 

First  we  need  a  well-known  tool  which  will  also  be  used  in  the  next 

section.  Let  g(0^,y),  i  *  0,  1,  be  functions  defined  over  a  subset  A 

of  the  real  numbers.  Assume  g(0 . ,y)  have  maximums  y.  in  A;  i.e., 

1  9i 


g(e*.yfl  )  “  “ax  g(04»y)  .  i  -  o,  1  . 

1  yeA  1 


Lemma:  If  for  every  y',y"  e  A,y“  >  y',  we  have 


gCej.y")  -  gCQj.y’)  <  g(®0,y")  -  gC0o,y’) 


then  <  y0Q' 


Proof:  We  have,  by  definition. 


gC^.y)  -  g(91,ye  )  <  o  ,  »yt  A  ,  i  -  0,  1. 


Suppose  y  >  y.  ,  y  «  A.  Then,  by  hypothesis,  with  y  -  y”,  y.  «  y' , 

0  eo 


g(®j.y)  -  g(VyeQ)  i  g<90»y)  -  g(eo,y0o^ 


<  o  . 


f 


g( 6! »y)  <  s^.y©  ) 


*  y  >  ye  • 

eo 


Therefore, 


y*  <  y 

♦l 


We  now  show 


Proposition  8;  y  ^  £  y^,  n  1. 

Proof:  Let  0  ■  n  +  1,  1  ■  0,  1,  so  that  in  the  lemma 

gC^.y)  •  log  (^(y)  -  p)  1  F(y)^^ 

-  log  (^(y)  -  n)  -  log  F(y)  +  log  (1  -  F  *H*i(y))  , 

i  *  0,  1  . 


Let  A  be  the  set  of 
are  assuming  v^^  >  vQ, 
tains  yn+i,  i  -  0,  1. 


y  for  which  ♦^(y)  -  p  >  0.  Since,  throughout,  we 
A  is  non-empty  and,  by  (7)  and  Proposition  7,  con- 
Hovever,  for  every  y"  >  y'(y',  y"  e  A)  we  have 


g^.y")  -  gCe^y’) 


log 


l-F°+1(y')‘ 


i  ■  0,  1 


For  (8)  to  hold  it  is  sufficient  that 


log 


l-F^y") 

l“Fn+1(y’) 


<  log 


i-rn(y") 

1-Fn(y’) 


*  y"  >  y’  , 


or  equivalently, 


1-F“(y’) 


1-Fn(y") 


This  inequality  will  hold  if 


1-Fn(y) 


is  non-increasing  in  y,  or  if 


1-Z 


n 


is  increasing  in  Z  for  0  <  Z  <  1.  However,  for  n  >  1, 


1-Z 


n+1 


n-1 


n  x.  /  n~ 

-  E  z  /  I 


1-Z  k-0  k-0 


/n-1  , 

E  z 

k-0 


-  1  ♦  (  I  -\)_1  . 

k-1  ZK 


which  is  increasing  in  Z. 

By  Proposition  8  lim  {y  }  exists.  We  denote  this  limit  by  y*. 
n-*»  n 

The  expression  <J^(y)  ~  U  represents  the  difference  in  the  expected 
life  if  one  spare  is  available  compared  to  having  none  available  when  the 
replacement  is  made  at  time  y.  In  our  notation,  y^  -  x^  is  the  optimal 
time  to  make  such  a  replacement.  We  have  as  a  result  of  Proposition  8  a 


Corollary:  ♦,  (y  )  “  w  is  non-increasing  in  n. 


•  Relationship  Between  Optimal  and  Optimal  Equal  Interval  Schedules 


Proposition  9; 


lim  <1*  “  lim  v 
n-Ke  n  n>®  n 


*  v 


Proof;  Consider  the  case  of  x*  >  0.  From  (1)  and  (2)  and 
Proposition  5  on  letting  n  ®,  we  have 


x* 

v  -  /  F(t)dt  +  F(x*)v  , 
0 


l*e«  | 


v 


f  F(t)dt 
6  F(x* )  * 


However,  from  (7) 


lim  (x)  -  p 
n 


{/  F(t)dt  -  F(x)p},  -v  x  . 


Thus,  in  particular. 


lim  <|>  >  lim  <j>  (x*) 
n  ■+«  n  n-*=°  n 


x* 

-  /  F(t)dt/F(x*) 
0 


On  the  other  hand,  since  v  >  <b  , 

—  n 


v 


Now  suppose  x*  -  0.  Ue  have  by  (7)  that 


lim  4>  >  f*  F(t)dt/F(x)  ,  v 
n->®  n  ~  a 


hence. 


x 

lim  4*  >  lim  /  F(t)dt/F(x) 
n-*»  n  —  x+0  o 


-  l/f(0)  . 


If  11m  v  ■  ®,  then,  by  Proposition  5,  f(0)  -  0. 

nx»  n 

11m  v  <  <*>,  by  Proposition  4, 

QXa  n 


1/v  *  lim  1/v 

n-x”  n- 1 


■  lim  r(x  ) 

n-x»  n 


“  r(0) 


-  f(0)  . 

Thus, 


lim  4<  >  v 
n-*»  n 


Since  lim  <b  <  v  »  v,  the  equality  follows  again 
n-x»  n  —  n 

Let 


<Ky)  -  lim  4»  (y) 


-  /  F(t)dt  +  F(y)v  ,  y  >  0 
0 


Since  4<n(y)  *  <Ky),  y  >  0,  and  <i>n(y)  and  <Ky)  are  continuous, 
then 


<Ky*)  *  max  (Ky) 
y>0 


if  y*  >  0.  If  v  <  ®,  we  have,  for  the  same  reasons,  that 


$(x*)  *  max  $(y) 
y>0 

“  v  , 


implying  that 


F(t)dt 

F(x*) 


v 


By  applying  the  lemma  of  Section  3,  with 


g(0o»  y)  “  log  (4<n(y)  -  V) 


and 


gC^.  y)  •  log  (<Ky)  -  n)  , 


in  the  same  manner  as  in  the  proof  of  Proposition  8  we  can  show  that 
y**  >  0  and  <Ky**)  ■  max  <Ky),  then  y**  c  A  and 


V  n 


Proposition  10: 


y*  • 


Proof:  If  x*  >  0,  y*  >  0,  then  by  the  preceding  statements  x*  and 
y*  are  each  the  largest  values  that  maximize  <Ky)  over  (0,<*>).  Thus, 
x*  *  y*.  We  need  only  show  that  X*  and  y*  are  positive  or  zero 
together.  Suppose  x*  >  0.  Then  let  y**  -  x*  in  the  remarks  preceding 
Proposition  10  implying  y*  y**  ■  x*  >  0.  Suppose  y*  >  0.  Then 

l/f(0)  <  <Ky*)  <  -  • 
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<Ky*)  >  1/f (0)  , 


then  x*  >  0  since,  by  Proposition  4,  x*  ■  0  implies  v  ■  l/f(0). 

<Ky*)  -  1/f (0)  , 

then  v  -  1/f ( 0) .  However  x*  is  the  largest  value  such  that 

(Kx*)  “  max  <t>( y ) 
y>0 

-  l/f CO)  . 


Therefore  x*  >  y*  >  0. 


Example : 

Uniform  (0,1). 


’  2  + 

-  *  (1  +  <1  -  (1  -  y)n)(l  -  y)>  . 


Thus,  y  is  that  value  of  y  that  maximizes 
’  n  J 


(1  -  y)(l  -  (1  -  y)n)  . 


If 


Letting  Z  ■  1  -  y,  we  see  that 


5.  The  Expected  Number  of  Spares  Ueed 

While  there  may  be  n  spares  available  to  extend  the  life  of  the 
system,  because  of  possible  failure  of  the  system  before  all  spares  are 
used,  the  actual  number  used  is  a  random  variable  depending  on  F  and  the 
schedule. 

Let  S  ■  (Z  ,  Z  ,  . ..,  Z, ,  T)  denote  the  schedule.  Let  u 
n  n-i  l  n 

denote  the  expected  number  of  spares  used.  Then,  u^  satisfies  the 
recursive  relationship 

(9)  u  -  F(Z  )(1  +  u  ,)  ,  n  >  1 

n  n  ni  — 

u0  -  °  . 

In  case  S  is  the  equal  interval  schedule  with  interval  length  y, 

then 


(10) 


un  -  F(y)  +  F(y)un_1 


_  F(y)-Fn~>‘1(y) 

F(y) 


> 


which  can  also  be  obtained  from  the  geometric  distribution. 


Examples 


(a)  Uniform  (0,1) 


If  the  optimal  schedule  is  used  then  (9)  yields 
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u 

n 


(1  -  x  )  +  (1  -  x  )  u 

n  n  n  i 


Vi'1 


+  Vi> 


If  the  optlnal  fixed  Interval  schedule  is  used,  then  (10)  yields 


u 

n 


f-Li1'"  .  rJLi 

^ir+l ' 


n+1  - 


(1-1/n+l) 


1/n 


n+1 

n 


l-(l/n+l)1/n 


r  1  U/n  n 
'■trU^  n+1  . 

l-(l/n+l)1/n 


These  values  will  be  tabulated  in  Table  1. 


(b)  Two  exponentials  in  parallel. 

In  this  case,  the  formula  (9)  for  the  optimal  schedule  yields, 


un  -  (1  -  (2v  -  1)“2)(1  +  u  ,) 

n  n  l  n  i 


These  values  will  be  tabulated  in  Table  2.  We  shall  be  especially 
Interested  in  comparing  these  values  with  expected  number  of  components 
used  in  the  situations  described  in  Sections  6  and  7. 


A 

.  J 

•  'J 


M 
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J^L 


6.  Systems  With  More  Than  One  Replaceable  Component 

In  the  previous  sections  it  was  assumed  that  when  a  replacement  was 
scheduled,  the  entire  system  was  replaced.  However,  example  (b)  arose  in 
the  context  of  the  system  consisting  of  a  two  independent  component 
module,  each  component  having  the  same  exponential  life  distribution. 

This  suggests  that  efficiencies  can  be  gained  by  exploiting  the  structure 
of  the  system  and  the  individual  life  distributions  of  the  components. 

We  concentrate,  here,  on  the  specific  structure  of  example  (b)  and 
consider  schedules  that  call  for  the  replacement  of  individual 
components.  A  number,  n,  of  spare  components  (as  distinguished  from  spare 
systems)  are  on  hand.  The  objective  is,  still,  to  attempt  to  maximize  the 
expected  life  of  the  system. 

To  put  the  problem  in  a  general  context,  a  working  system  can  be 
thought  of  as  being  in  state  (s,t),  meaning  that  s  time  units  have 
elapsed  since  it  was  known,  either  by  replacement  or  deduction,  that 
component  1  was  in  working  order.  The  same  meaning  is  given  to  t  for 
component  2.  The  deduction  that  component  1  is  working  is  made  when  the 
inspection  (which  we  assume  takes  place)  of  component  2  after  it  is 
replaced  reveals  that  it  has  failed,  implying  that  component  1  is  working, 
and,  because  of  its  exponential  character,  is  as  good  as  a  new  component. 
In  state  (s,t),  at  least  one  of  the  components  is  working. 

Remark;  For  the  above  to  be  applicable  we  are  assuming  that: 

a)  inspection  of  a  component  is  destructive  (for  otherwise,  one  would 


continue  to  use  a  working  component  rather  than  replace  it). 


b)  inspection  is  informative  in  the  sense  that  it  tells  us  whether  or 


not  the  component  was  functioning  immediately  before  Inspection. 

Let  w  (s,t)  denote  the  expected  life  of  the  system  when,  Initially, 
the  state  is  (s,t),  n  spares  are  available,  and  an  optimal  replacement 
schedule  is  used. 


Special  Case:  Assume  (s,t)  ■  (0,0)  and  n  *  1,  and  assume  that 
component  1  is  optimally  replaced  at  time  Z^.  Then 


w^(0,0) 


max  {2  /  ye  ^(l-e  ^)dy  +  e  2Z( Z  +  3/2)  +  (l-e  Z)e  Z(Z  +  3/2) 
z  0 

+  e"Z(l  -  e_Z)(z  +  1)} 

Z 

/ 

0 


-Z 1  -z, 


-2  /  1  ye  y(l  -  e  y)dy  +  (5/2  +  Z  )e  1  -  (l  +  Z  )e 

A  * 


-2Z1 


It  is  seen,  in  a  straight  forward  manner,  that 


_Z1 

e  1  -  1/2 


or 


Zx  -  log  2  , 


from  which  it  follows  that 


Wl(0,  0)  -  13/8  . 


The  Case  of  n  >  1. 


When  n  >  1,  the  problem  of  obtaining  an  optimal  component  replace¬ 
ment  schedule  appears  to  be  difficult  to  solve.  For  the  purpose  of,  at 


least,  achieving  lower  bounds  on  what  could  be  accomplished  by  an  optimal 
schedule,  we  restrict  ourselves  to  a  subclass  of  schedules  for  which  the 
analysis  is  feasible. 

We  assume,  initially,  that  (s,t)  -  (0,0)  and  consider  schedules  of 
the  following  form:  when  n  2,  at  time  xq,  component  1  is  replaced  and 
then  inspected.  If  the  inspection  reveals  that  the  component  was  working, 
then  component  2  is  also  replaced.  If  the  Inspection  shows  that  the 
component  was  not  working,  then  component  2  is  left  intact.  When  n  -  1, 
component  1  is  replaced  at  time  . 

The  rationale  for  the  type  of  rule  is  based  on  the  exponential  life 
distribution  of  the  components,  the  parallel  structure  of  the  system,  and 
mathematical  convenience.  If,  when  the  system  is  working,  one  component 
is  not  functioning,  then  the  other  must  be,  and,  since  it  is  working  it 
must  be,  statistically  speaking,  equivalent  to  a  new  component.  Both  com¬ 
ponents  are  replaced  when  the  one  inspected  is  functioning  so  that  the 
state  is  always  (0,0)  after  a  scheduled  replacement.  However,  the 
double  replacement  is  not  completely  inefficient,  for  if  the  one  inspected 
is  working  there  is  an  enhanced  possibility  (depending,  of  course,  on 

x  )  that  the  other  will  have  failed. 
a 

Let  vR  denote  the  maximal  life  of  the  system  when  n  spare  com¬ 
ponents  are  available  and  the  class  of  schedules  is  restricted  to  the 
above  stated  class.  Then  v  satisfies 


n-1  n-2 


If  n-1,  then  x^  -  x^  -  log 


2. 


On  Substituting  x  In  (11)  we  get,  after  some  calculations,  that 

U 


n  2 


1{3  +  (Vl+Vn-2"2) 


2v  -1 
n-1 


■}  ,  n  >  2  . 


To  calculate  un>  the  expected  number  of  spare  components  used  In  the 
life  of  the  system,  we  have 


u0-  ° 


-XU2 


-  1  "  (1  -  •  ) 

-  3/A 


X  X 

un  ■  “  e  **)•  0(1  +  un-i>  +  e  11(2  +  Un-2J  »  n  >  2 


-V 


The  formulas  for  x  ,  v  ,  and  u  are  recursive  and  are  readily  cal- 

n  n  n 

culated  numerically.  (See  Table  3.)  In  comparing  x  ,  v  and  u  with 

n  n  n 

**  **  A* 

x  ,  v  and  u  It  should  be  kept  In  mind  that  n,  in  the  former  case, 
n  n  n 

because  the  entire  system  is  replaced  is,  in  fact,  twice  the  number  of 
components  of  the  latter  case. 


.  The  Option  of  Non-destructive  Inspections 


In  Section  6  it  was  assumed,  presumably  due  to  destructive 
inspections,  that  an  inspection  implied  a  replacement.  Here  we  assume 
that  Inspections  are  not  destructive  so  that  if  a  component  is  found  to  be 
working,  it  is  continued  in  use.  Without  assuming  inspection  costs, 
continuous  surveillance  would  be  optimal  leading  to  an  expected  system 
life  equal  to  — y  +  1.  Instead  of  introducing  inspection  costs  into  the 

model  we  shall  derive  formulas  for  computing  the  values  of  the  charac¬ 
teristics  assuming  the  schedule  derived  in  Section  6.  In  this  way  we  can 
compare  destructive  and  non-destructive  Inspections  with  respect  to  a 
fixed  schedule. 

Let  u*  denote  the  expected  number  of  replacements  made  during  the 
n 

life  of  the  system.  Then 


i*€t  | 


u* 


X. 

(1  -  e  l) 


0  +  e 


-2xl 


U* 


x, 

2e  1{l  - 


-xli 


-X1 


-xl. 


„  2e  (1-e  ») 


1-e 


-2xl 


-  2/3  . 


For  n  >  2, 


u*  -  (l  -  e  n)2 


0  +  e 


u*  +  2e"  n(l  - 

fl  v 


l)(l 


+• 


Let  v*  denote  the  expected  life  of  the  system.  Then 
n 


V.  -  3/2 


2  I*"  +(!-(!-  e-Xn)2)  I 


LA  A  A 

+  e“  n  v*  +  2e~  n(l  -  e“  n)  v*  .  ,  n  > 

n  '  n-1  — 


i.e.,  for  n  >  1, 


i-S-  I  f”  Cl  -  (1  "  «"V)dt  +  2e“Xn(l  -  e_Xn)v*n-1} 


2x  0 
—  n 


i _  (3  .  2  -*n  e~  n  +  2e'  n(l 

~  *2  e  2 


-  •'  ">  vS-il 


Values  of  v*  and  n*  can  be  recursively  calculated  and  are  tabulated  In 
n  n 

Table  3. 

Considering  that  n  In  Table  2  Is,  In  fact,  twice  the  value  of  n 
In  Table  3,  replacement  of  components  rather  than  systems  produces  a 
substantial  gain  in  expected  life  for  a  given  number  of  spare  components. 
On  the  other  hand,  in  comparing  the  two  parts  of  Table  3,  non-destructive 
rather  than  destructive  inspection  produces  only  a  modest  gain. 


»*****-•- 
■s.'  V  s*  • 


■  «.  *  .*•  m'*  W***  »**  «  *  t  *  •  .  * 

**a  *  .  *  "  4 '  *  • 


s’ 


8.  Asymptotic  Distribution  of  System  Life 

Let  Pn(t)  denote  the  probability  that  system  life  exceeds  t  vtoen 

n  spares  are  available  and  an  optimal  replacement  schedule  is  used.  Let 

P  (t)  denote  the  same  probability  vftien  an  optimal  equal  interval  replace- 
n 

ment  schedule  is  used.  When  x*  *  y*  ■  0  and  v  ■  <J>  »  l/f(0)  >  0,  it 
would  seem  that  the  limiting  distribution  as  n  -*■  ®  of  system  life  should 
be  exponential  with  rate  f(0).  That  this  is  true  depends  on  the 
following. 

2 

Proposition  11; 


l  x,  -  lim  ny 
j.  ^  i  rv+“  n 


Proof : 

true.  Also, 


In  the  cases  where 
since 


x* 


y*  >  0,  the  proposition  is  trivially 


n 

v  <  l  x  + 
H-J.l  i 


and 


<l»  <  ny  +  p  , 

n  —  n 


the  proposition  is  immediate  when  v  ■  Thus,  consider  the 


remaining  case  vftiere  v  •  ■  l/f(0)  <  «°  and  x*  ■  y*  *  0.  Since  the 


^e  are  indebted  to  D.  R.  Smith  for  some  remarks  which  were  very  helpful 
In  establishing  this  proposition. 


sequences  {x  }  and  {y  }  are  non-increasing  and  because  f  is 
n  n 

continuous  of  x  ■  0,  for  any  e  >  0  there  exists  an  N  such  that  for 
n  >  N 

(12)  F(x  )  <  (f(0)  +  e)x  ,  F(y  )  <  (f(0)  +  e) y  . 

n  -  n  n  -  n 

On  Iterating  (5) 


n-1 

v  -  v  .  >  II  F(x  )(v  -v.)  ,  *  n  2  . 

n  n-1  J  10 


Since  v  <  ®,  lim  (v  -v  , )  ■  0  implying 
A  -+ao  n  n- 1 


n  f(x  )  -  o  , 
J-i  J 


which,  with  (12),  implies  x^  *  «.  From  (7),  with  y  -  yn,  on 

letting  n  -*■  <*>  one  has  that 


lim  Fn(y  )  “  0  , 

n+®  “ 


which,  with  (12),  Implies  lim  ny  •  ®  . 

n-Ko  n 

We  now  have 


Proposition  12;  If  x*  -  0,  then  for  every  t  >  0, 

lim  P  (t)  -  lim  P  (t) 
n-*=  n  n+*°  n 

-  e-f(°)t 


if  f(0)  >  0 
if  f(0)  -  0 


Hence,  for  every 


nOO 

Proof:  Fix  t.  From  Proposition  11,  xj  “  °°* 

k  there  exists  an  N  such  that  for  n  >  N,  J Xj  >  t.  Thus,  since 

the  x  *8  are  decreasing  the  component  in  use  at  time  t  or  before  is 

n 

never  older  than  x^.  Denote  by  r^(t),  the  failure  rate  function  for  the 
system  life.  Then 

min  r(s)  <  r  (t)  <  max  r(s)  ,  0  <  t  <  t  . 

0<s<xk  —  n  —  0<s<Xk  ~  ~ 


Thus,  using 


Pn(t) 


JZr  (x)dt 
-0  n 


-t  max  r(s)  -t  min  r(s) 

0<s<x,  (Kb<x, 

e  k  <  Pn(t)  <e  k  ,  v  k  and  n  >  N 


Letting  n  -*  ®  and  then  k  -*■  ®  and  since  x^  0,  we  have  the  value  for 
lim  P  (t).  The  value  for  lim  P  (t)  is  proved  similarly. 


9.  A  Multi-Failure  Generalization 

Let  us  now  suppose  that  it  takes  1  (i  1)  vital  component  failures 
to  cause  a  system  failure.  Component  replacement  occurs  at  a  component 
failure  or  when  scheduled.  In  previous  sections  i  *  1.  As  before  we 
suppose  there  are  n,  n  i-1,  spares  and  we  let  v(n,i)  denote  the 
expected  life  of  the  system  under  an  optimal  schedule  when  the  system 
fails  at  the  moment  of  the  i**1  vital  component  failure  and  n  spares 
are  available.  Then  v(n,i)  satisfies  the  optimality  equation 


v(n,i) 


max 

0<x<T 


x 


{/ 

0 


F(t)dt  +  F(x)v(n-l,i)  +  F(x)v(n-l,i-l)  }  ,n  i,i  1, 


v(i-l,i)  *  ip.  (Replacement  only  occurs  at  component  failures  when 
n  -  i-1.)  , 

v(n,  0)  -  0  ,  n  ^  0  . 

Let  us  denote  by  x(n,  i)  the  largest  value  that  maximizes  the  above  ex¬ 
pression  within  brackets.  The  following  proposition  whose  proof  by  induc¬ 
tion  on  k  =  i  +  n  will  not  be  given  as  it  is  almost  Identical  with  the 
proof  of  Lemma  3  of  [ 1 ] . 


Proposition  13: 

(a)  v(n  +1,  i  +  1)  -  v(n  +1,  i)  >.  v(n,  i  +  1)  -  v(n,  i) 

(b)  v(n  +  1,  i)  -  v(n,  i)  >  v(n  +  2,  i)  -  v(n  +  1,  i) 

(c)  v(n,i+l)  -  v(n,i)  v(n,i+2)  -  v(n,i+l) 

From  the  lemma  of  Section  3  and  the  above  proposition  we  have  the 


following. 


10.  Further  Remarks 


Proposition  3  and  13(b)  are  relevant  to  the  question  of  how  many 

spare  parts  should  be  stocked.  The  propositions  assert  that  the 

incremental  gain  in  expected  life  is  decreasing  in  n.  Thus,  if  the 

reward  associated  with  a  unit  of  expected  life  can  be  compared  to  the  cost 

of  a  spare,  there  is  a  determinable  N  such  that  it  is  not  economically 

efficient  to  shock  n  >  N  spares.  Although  it  seems  a  reasonable 

conjecture,  we  have  not  been  able  to  show  that  <l>  -  <p  ,  i  in  n. 

n  n-1 

Suppose,  instead  of  maximizing  expected  life,  it  is  more  relevant  to 
maximize  the  probability  that  the  system  survives  t  units  of  time. 

Then,  in  the  case  of  i  ■  1,  if  the  intervals  between  scheduled  replace¬ 
ments  of  n  spares  is  (z  ,  z  . . z, } ,  then 

n  n-1  1 

P(L  >  t)  »  5  F(z.)  , 

i“0  * 

where 

n 

(13)  l  z  »  t  ,  z  >  0,  i  -  0,  ...,  n  , 

i-0  i  1 

presuming  it  is  always  optimal  to  schedule  the  use  of  all  spares. 
Maximizing  P{L  >  t}  subject  to  (13)  is  equivalent  to  maximizing 

log  F(z^)  subject  to  (13),  a  classical  allocation  of  resources 

problem.  If  F  is  IFR,  then  log  F(z)  is  concave  and  it  is  well-known 
that  the  optimal  values  of  z^,  ...,  eq  are  equal;  i.e.,  z^  "  t/n+1, 
k  -  1,  ...,  n.  If  F  is  DFR  then  log  F(z)  is  convex  and  an  optimal 
schedule  is  to  have  z  »  t,  z,_  -  0,  k  ■  1,  . . . ,  n-1. 


TABLE  1.  UNIFORM  (0,1)  DISTRIBUTION 


n 

V 

X 

u 

y 

0 

.5000 

.0000 

.5000 

1 

.6250 

.5000 

.5000 

.6250 

.5000 

2 

.6953 

.3750 

.9375 

.6924 

.4226 

3 

.7417 

.3046 

1.3471 

.7362 

.3700 

1 

4 

.7750 

.2582 

1.7409 

.7674 

.3312 

1 

5 

.8003 

.2249 

2.1244 

.7911 

.3011 

1 

6 

.8203 

.1996 

2.5007 

.8098 

.2769 

2 

7 

.8364 

.1796 

2.8716 

.8250 

.2570 

2 

8 

.8498 

.1635 

3.2384 

.8377 

.2401 

2 

9 

.8610 

.1501 

3.6019 

.8484 

.2257 

3 

10 

.8707 

.1389 

3.9627 

.8576 

.2132 

3 

11 

.8790 

.1292 

4.3212 

.8656 

.2022 

3 

12 

.8864 

.1209 

4.6779 

.8727 

.1924 

3 

13 

.8928 

.1135 

5.0329 

.8789 

.1837 

4 

14 

.8985 

.1071 

5.3865 

.8845 

.1758 

4 

15 

.9037 

.1014 

5.7389 

.8896 

.1687 

4 

16 

.9083 

.0962 

6.0902 

.8942 

.1622 

4 

17 

.9125 

.0916 

6.4405 

.8983 

.1563 

5 

18 

.9163 

.0874 

6.7900 

.9022 

.1509 

5 

19 

.9198 

.0836 

7.1387 

.9057 

.1458 

5, 

20 

.9230 

.0801 

7.4867 

.9089 

.1412 

5, 

30 

.9449 

.0566 

10.9384 

.9315 

.1081 

7, 

40 

.9570 

.0438 

14.3581 

.9445 

.0886 

10, 

50 

.9647 

.0358 

17.7597 

.9531 

.0756 

11, 

60 

.9701 

.0303 

21.1494 

.9592 

.0662 

13, 

70 

.9740 

.0262 

24.5309 

.9638 

.0590 

15, 

80 

.9770 

.0231 

27.9062 

.9674 

.0534 

17. 

90 

.9794 

.0207 

31.2767 

.9703 

.0488 

19. 

100 

.9813 

.0187 

34.6435 

.9727 

.0451 

20. 

.0000 

.5000 

.9106 

.2768 

.6150 

.9336 

.2374 

.5296 

.8122 

.0869 

.3548 

.6167 

.8734 

.1255 

.3734 

.6176 

.8583 

.0959 

.3306 

.5627 

.7922 

.9797 

.0263 

.9836 

.8699 

.7023 

.4907 

.2422 

.9621 
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TABLE  2 


f(t)  -  2e“t(l-et) 


n 


0 

1.5000 

1 

1.7500 

2 

1.9500 

3 

2.1224 

4 

2.2765 

5 

2.4172 

6 

2.5476 

7 

2.6697 

8 

2.7849 

9 

2.8943 

10 

2.9987 

11 

3.0988 

12 

3.1950 

13 

3.2877 

14 

3.3774 

15 

3.4643 

16 

3.5486 

17 

3.6306 

18 

3.7105 

19 

3.7883 

20 

3.8644 

30 

4.5454 

40 

5.1259 

50 

5.6406 

60 

6.1078 

70 

6.5387 

80 

6.9407 

90 

7.3188 

100 

7.6770 

-  .0000 

1.0986  .7500 

.8472  1.4700 

.7191  2.1763 

.6370  2.8746 

.5785  3.5676 

.5339  4.2570 

.4984  4.9435 

.4693  5.6279 

.4448  6.3105 

.4238  6.9917 

.4056  7.6717 

.3896  8.3507 

.3754  9.0289 

.3626  9.7063 

.3510  10.3830 

.3405  11.0592 

.3310  11.7348 

.3221  12.4100 

.3140  13.0847 

.3064  13.7591 

.2524  20.4878 

.2196  27.1994 

.1970  33.9015 

.1802  40.5972 

.1671  47.2887 

.1565  53.9769 

.1477  60.6626 

.1402  67.3464 
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TABLE  3 


.6931 

.6931 

.6000 

.5083 

.4644 

.4195 

.3921 

.3647 

.3454 

.3266 

.3122 

.2983 

.2869 

.2761 

.2669 

.2582 

.2505 

.2433 

.2368 

.2307 

.1878 

.1623 

.1449 


1.5000 

.0000 

1.5000 

1.6250 

.7500 

1.8333 

.6666 

1.7812 

1.4375 

1.9166 

1.1666 

1.8858 

2.1127 

2.0536 

1.7273 

2.0013 

2.8136 

2.1654 

2.3381 

2.0930 

3.4753 

2.2729 

2.9437 

2.1883 

4.1723 

2.3670 

3.5500 

2.2707 

4.8328 

2.4583 

4.1710 

2.3537 

5.5237 

2.5415 

4.7807 

2.4289 

6.1858 

2.6222 

5.4079 

2.5035 

6.8714 

2.6975 

6.0224 

2.5730 

7.5353 

2.7707 

6.6524 

2.6415 

8.2169 

2.8401 

7.2715 

2.7062 

8.8822 

2.9075 

7.9033 

2.7699 

9.5607 

2.9721 

8.5263 

2.8309 

10.2272 

3.0351 

9.1594 

2.8907 

10.9033 

3.0958 

9.7856 

2.9483 

11.5706 

3.1550 

10.4197 

3.0049 

12.2450 

3.2124 

11.0486 

3.0597 

12.9128 

3.2685 

11.6837 

3.1135 

13.5859 

3.3231 

12.3149 

3.5950 

20.2823 

3.8116 

18.6817 

4.0056 

26.9702 

4.2266 

25.0900 

4.3699 

33.6533 

4.5938 

31.5261 

4.7005 

40.3332 

4.9267 

37.9823 

5.0055 

47.0109 

5.2334 

44 . 4542 

5.2900 

53.6870 

5.5193 

50.9386 

5.5577 

60.3619 

5.7881 

57.4334 

5.8112 

67.0359 

6.0426 

63.9370 

Acknowledgement : 

Ve  are  indebted  to  A.J.W.  Derman  for  the  micro-computer  programming 
required  to  compile  Tables  1,  2  and  3. 
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